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by Charles L. Zola and Laurence H. Flshbach 
Lewis Research Center 

SUMMARY 

Propellant consumption data are presented for constant-thrust, low-acceleration, 
Earth-Mars round-trip trajectories including escape and capture spirals at both planets. 
Calculations are based on a previous analysis in which the calculus of variations was ap- 
plied to the problem of minimizing the propellant consumption of the round-trip trajectory 
treated as a single unit. Solutions given are for initial accelerations between 1. 0 and 

O 

4. 0x10 meter per second squared and include optimum coast phases. Results for short 
trip times between 340 and 460 days and a long trip time of 1000 days are presented with 
various wait times at Mars. Also included are examples of mission profiles where the 
first and/or last Earth spirals are deleted. A simple example is given showing how data 
of the type presented may be used for mission analysis. The mission profiles and their 
data are limited in scope because of computational difficulties. These results, there- 
fore, are not sufficient for a complete study of the Mars round-trip mission. 


INTRODUCTION 

Electric propulsion will be attractive when it can provide lower vehicle gross weight 
for a space mission than can be achieved by an alternative propulsion method. Although 
the specific impulse of an electric rocket is high when compared with that of a chemical 
or nuclear rocket, propellant may still constitute a significant fraction of a given vehicle. 
It is thus important that the thrust vector during the long periods of propulsion be so di- 
rected as to minimize propellant requirements. This may be done by applying the calcu- 
lus of variations or other optimization theory in the calculation of the trajectory. 

Optimization theory has been used extensively in producing trajectory data for low- 



thrust one-way trips to various planets. Some analyses have assumed a vehicle with con- 
stant power but variable thrust (refs. 1 and 2) in which the thrust vector is optimally con- 
trolled in both magnitude and direction. Others interested in the performance of constant- 
power electric vehicles with constant thrust, typical of state-of-the-art engines, have 
formulated calculus of variations solutions for these and have given examples of such tra- 
jectory solutions over a range of travel times (ref. 3). For the constant-thrust problem, 
the thrust and total mass flow rate are held constant at prescribed values, while the varia- 
tional technique determines the optimum thrust direction programing and the placement 
of a coasting period (if any) to result in minimum propellant consumption. 

There has been considerable interest in analyzing the use of electric propulsion in 
round-trip interplanetary missions. Reference 4 gives an example of the use of constant 
thrust but does not use optimum control of the thrust vector. Another study (ref. 5) uses 
the variational calculus, but for the case of variable thrust. It is the purpose of this 
report to present examples of short- and long-duration round trips to Mars for the 
constant -thrust vehicle with the thrust vector direction optimally controlled. 

Many pairs of minimum propellant outbound and inbound trajectories would suffice to 
perform a given round trip, but total propellant requirements for the overall mission 
could be overestimated unless the best trajectory pair is found. The best pair was found 
in reference 5 by a direct trial and error method that combined precalculated one-way 
trajectory data into optimum round trips. In the case of constant thrust, however, the 
additional parameters, thrust and total mass flow rate, would require an almost prohibi- 
tive amount of precalculated data to employ a trial and error method. Therefore, the 
data presented in this report are based on a different approach involving a further appli- 
cation of the calculus of variations. The basic technique used is fully described in refer- 
ence 6 where the theory developed in reference 3 for one-way trajectories is extended to 
the problem of optimizing the complete low-thrust round trip. Reference 6 shows that 
certain parameters of the variational solution relate the optimum pair of outbound and 
return heliocentric trajectories. The implicit relations derived in this reference are 
sufficient to allow the appropriate pair of trajectories to be readily identified and calcu- 
lated simultaneously. 

The basic computer code of reference 3 and the analytic method of reference 6 (re- 
duced to two dimensions) were combined to provide the trajectory data of this report. An 
idealized circular orbit planetary model is assumed. Each round-trip trajectory is cal- 
culated straightforward from beginning to end. Every converged solution represents a 
minimum propellant round trip for a given thrust level and total mass flow rate. Mission 
time and waiting time are, however, dependent variables. External controls must be 
used to keep mission time and wait time constant when other parameters are varied. 

Although this procedure yields accurate trajectory data, within the confines of the 
model, it has been found in preparing the limited amount of data for this report that cal- 
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Figure 1. - Nominal mission profile for manned Mars mission using 
electric propulsion. 

culation difficulties for certain cases can result in excessively long computation time, in 
many instances 20 minutes or greater on an IBM 7090. For this reason it becomes pro- 
hibitively expensive to run enough trajectory data to cover all the cases necessary for a 
thorough mission analysis. Nevertheless, the data presented here can be utilized in pre- 
forming examples of mission analyses employing electric propulsion (appendix B). 

The chief value of the data lies in their utilization as a check on more approximate 
but faster trajectory schemes and as an illustration of the characteristic behavior that 
can be expected. For example, the trajectory method used in the mission analysis of 
reference 7 was initially checked by this method. 

Minimum-propellant, constant-thrust, Earth-Mars round-trip trajectories are pre- 
sented in the low range of trip times from 340 to 460 days with various wait times in a 
1. 1 Mars radii orbit. Data are also included for a trip time of 1000 days with waiting 
times of 300, 400, and 500 days. 

It is shown in appendix B that the effects on payload of unionized propellant leaving 
the vehicle and inefficiencies associated with ionizing the propellant can be determined 
without repeating the trajectory calculations. 


ANALYSIS 

A typical round-trip trajectory profile to Mars is shown in figure 1. The problem 
model assumes a two-dimensional solar system with the planets in circular orbits about 
the Sun. The round-trip trajectory is approximated as a "patched” sequence of two-body 
planetocentric and heliocentric phases. The assumed constants of the planetary model 
are the following: 
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Earth 


Mars 


Distance from Sun, m 

Radius of planet, m 

Gravitational constant, m /sec . . . . 

Parking orbit in planet radii 

Angular velocity about Sun, radians/sec 


. . 1.495X10 11 

. . 6.37123X10 6 
3. 99 1666 IX 10 14 

1.1 

1. 9909836x10 ^ 


2.2779179X10 11 
3. 33215X10 6 
0. 4297866 lxlO 14 
1 . 1 

1. 0585766x10 ^ 


Heliocentric Trajectories 

The heliocentric trajectories (1-2 and 5-6) in figure 1 begin and end with the vehicle 
in circular orbit about the San at each respective planet- Sun radius. During heliocentric 
flight, the optimal thrust vectoring program and placement of a coasting phase is governed 
by the variational calculus. Each trajectory is solved by using the basic computer code 
discussed in reference 3. The code incorporates a Runge-Kutta numerical integration 
technique with step size control to limit truncation error. The integration procedure is 
coupled with a three-variable Newton-Raphson iteration scheme to solve the boundary 
value problem presented by each heliocentric trajectory. 


Planetocentric Trajectories 

Planetocentric escape and capture phases (0-1, 2-3, 4-5, and 6-7) are included in 
the model because their effect on the mission as a whole can be important. Two means 
of executing escape and capture maneuvers are included in the model, and in this report, 
because their effect on the mission as a whole can be important. 

The first escape method is a continuous low-thrust spiral between circular orbit at 
1. 1 planet radii and escape energy relative to the planet. Escape is assumed to occur 
sufficiently far from the planet that the vehicle can be considered out of the planet's grav- 
ity field and that it has zero velocity relative to the planet. At this point heliocentric 
flight begins. With long propulsion times, low-thrust planetocentric maneuvers should 
use optimum or close to optimum thrust vector control. Reference 8 shows that constant 
tangential thrust acceleration (parallel to the velocity vector) agrees within 1 percent of 
the variational solution. The results for constant tangential thrust are expected to like- 
wise agree and are much easier to calculate. Even the simple tangential thrust spiral 
calculation would be time consuming; therefore, data from precalculated generalized tan- 
gential spiral solutions (ref. 4) have been fitted within 1-percent accuracy by empirical 
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curves. In this way, time and propellant requirements for escape or capture spirals can 
be found for any given pair of thrust and total mass flow rate. 

The second maneuver used a parabolic conic section to replace the low-thrust spirals 
necessary when electric thrustors are employed for planetocentric escape or capture. To 
escape, it is assumed that a high-thrust booster places the vehicle on an escape parabola. 
No penalty for performing this high-thrust maneuver has been included in the data of this 
report. The vehicle coasts until it can be assumed that it is in heliocentric space having 
zero velocity relative to the planet. The time spent coasting on this parabola has not been 
included in the mission time. At this point the electric thrustors are turned on and the 
heliocentric transfer begins. The use of parabolic capture trajectories represents high- 
thrust braking into elliptic or circular orbits or atmospheric braking followed by descent 
to the surface of the planet. 


Calculation of Complete Trip 


Combining the results of the planetocentric maneuvers with the calculus of variations 
solution to the heliocentric transfer yields the complete vehicle performance. The mass 
at any time M, the thrust F, and the total mass flow rate M T qt are suf fi c i ent to cal- 
culate both the trajectory and accompanying mass change. 

The M T Qrp is assumed to be inclusive of all mass flow rates regardless of whether 
or not they contribute to thrust. In presenting the results of this report, it has been 
found to be advantageous to use alternate parameters that are defined in terms of M T q T 
and F. These parameters are the thrust to initial mass ratio (initial acceleration) F/Mq 
and the effective power to initial mass ratio P/Mq. The F/Mq ratio is defined by the 
thrust and initial mass, and P/Mq can be expressed as 


P _ 1 /_ f \ F 
M 0 2 \ M 0/ ^TOT 


( 1 ) 


where F/M TOT is the effective exhaust velocity of the expelled mass. The P/M Q is 
shown in appendix B to be a useful parameter in that it is of the order of the actual beam 
power and total propulsive power required by the vehicle and differs only by thrustor effi- 
ciency functions. 

When P/Mq and F/Mq are used, the calculation of a typical round-trip trajectory 
commences with arbitrarily choosing an initial mass Mq and thus defining F and 
MpoT' Tlie M 0 * s identical to M j when the initial Earth spiral is replaced by a high- 
thrust maneuver. The mass change evaluated for each phase then determines the initial 
mass of each succeeding phase. The process is repeated throughout the trajectory and 
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Figure 2. - Final mass ratio against thrust to initial mass ratio. Mission time, 
380 days,- wait time, 10 days; four spirals. 


finally results in a terminal mass fraction referred to as Mp/Mg. The final mass 
corresponds to Mg for no terminal spiral or to when a terminal spiral is included. 

The trajectory calculations presented here do not include arbitrary (i. e. , nonpropel- 
lant) mass changes at Mars. In any actual round-trip mission, some arbitrary amount of 
mass (e. g. , the landing and exploration system) might be jettisoned or left behind at the 
start of the return trip. Final mass fraction of the trip Mp/Mg decreases as the frac- 
tion of jettisoned mass to initial spacecraft mass increases. 


RESULTS AND DISCUSSION 

Trajectory data are presented in the form of graphs (hereafter called performance 
maps), such as figure 2, showing final mass ratio Mp/Mg as a function of thrust to ini- 
tial mass ratio F/Mg. At each point on a map an effective specific impulse (effective 
jet velocity divided by the gravitational constant of Earth, 9. 80665 m/sec^) is defined 
since the power to mass ratio is directly proportional to the thrust to mass ratio and jet 
velocity (eq. (1)). 

Performance maps are presented for three classes of missions: short mission 
times with four spirals, short mission times with the initial and/or final Earth spirals 
omitted, and a long, four- spiral mission time with various waiting times at Mars. 

All the performance maps presented in this report (appendix C) have similar charac- 
teristics. The purpose of this section is to examine general changes in performance 
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brought about by changing the mission parameters, such as mission time and waiting 
time, and changing the profile by the omission of spirals. 

Figure 2 is typical of the performance maps, in this case a four- spiral trip for a 
380-day total mission time with a 10-day wait at Mars. Along each constant P/Mq line 
there exists an optimum F/Mq, which yields a maximum M^/Mq. The M^/Mq contin- 
ues to improve in the direction of higher P/Mq. For finite F/Mq, as P/Mq becomes 
infinite (infinite specific impulse), the mass ratio approaches unity. 

In figure 2, the circled point at the extreme left of each P/Mq curve is where the 
round-trip trajectory degenerates into all propulsion on both legs of the trip. For each 
P/Mq this point defines the lower bound on F/Mq. All-propulsion, or zero coast time, 
solutions result when the energy requirements are so severe that the thrust or must oper- 
ate continuously over the available time. To the right of the all-propulsion point at each 
P/Mq some coasting occurs on both outbound and inbound legs. 

Examination of a typical P/Mq line shows that as F/Mq is decreased (heading 
towards the all-propulsion point), M^/Mq starts to decrease rapidly. As was pointed 
out in the INTRODUCTION, computational difficulties arise when certain cases are run; 
typical of these are those near the all-propulsion boundary. Therefore, to save computa- 
tion time, solutions were usually stopped before the all-propulsion points were reached. 
The dashed line in figure 2 represents the F/Mq where the solution would have normally 
been terminated. 

The sketch shows that as F/Mq is decreased, both outbound and return coast times 
approach zero. The outbound coast time tends to level off until the return coast time is 
short, then both appear to approach zero simultaneously. This is the area of computa- 
tional difficulty mentioned previously. Solutions were stopped at an F/Mq that corre- 
sponds to about a 5-day outbound coast (indicated by dashed line). 

The bound that appears on the left-hand side of the performance maps (appendix C) 
represents the F/Mq where solutions were stopped and the data should not be extrapo- 
lated to lower F/Mq. 

Figure Cl(a) in appendix C presents a plot of a round-trip trajectory for a 30-day 
wait, 380-day trip with no initial or final Earth spirals. Figure Cl(b) presents a similar 

plot for a 400-day wait, 1000-day trip with four spirals. Ta- 
ble Cl in appendix C gives trajectory variables as functions 
of time for the 380-day trip. The purpose of presenting these 
data in this manner is to allow the reader to check any re- 
sults he has obtained with his own trajectory program against 
those presented here. 

The remainder of this report will be concerned with the 
effects of trip time, wait time, and mission profile on final 
mass ratio. These effects will be shown by plotting maximum 
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Figure 3. - Effect of mission time on maximum final mass Figure 4. - Effect of wait time on maximum final mass ratio. Four- 

ratio. Four-spiral missions; wait time, 10 days. spiral missions of 380 days. 


Mp/MQ against P/Mq as determined from the appropriate performance maps found in 
appendix C. 


Four Spiral Trips 

Figure 3 shows the effect of trip time on final mass fraction for trips with a constant 
wait time. The performance maps for trip times of 340, 380, 420, and 460 days with a 
10-day wait are shown in figure C2 of appendix C. As trip time is increased from 340 to 
460 days, the final mass fraction can be seen to increase from 0. 441 to 0. 634 at a P/Mq 
of 150 watts per kilogram. Longer trip time allows lower acceleration and, hence, 
higher specific impulse at the same power, which results in increased M^/Mq. 

Figure 4 shows the effect of wait times of 10, 30, and 50 days on final mass fraction 
for a 380-day trip. Figure C3 shows the performance maps for these missions. It can 
be seen in figure 4 that M f /Mq decreases with increasing wait time dropping from 
0. 511 at a 10-day wait to 0. 417 at a 50-day wait for P/M Q of 150 watts per kilogram. 
For a fixed mission time and an increasing wait time, this decrease in M f /Mq is a re- 
sult of forcing the transfer time to be shorter, which results in an effectively more diffi- 
cult mission. The optimum wait time is therefore zero days, but the choice of wait time 
is not arbitrary and will be dictated by the time required to complete the objectives of the 
mission. 
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Figure 5. - Effect of omission of spirals at Earth on maxi- 
mum final mass ratio. Mission time, 380 days; wait 
time, 30 days. 


Omission of First and/or Last Spirals 

Maps are presented in appendix C (figs. 

C4 and C5) for 340- and 380-day trip times, 
respectively, with a 30-day wait in which one 
or two spirals have been omitted. Omission 
of the first spiral corresponds to a boost of the 
electric vehicle to escape energy through an- 
other type of rocket. Figure 5 shows the effect 
on final mass fraction of removal of the first 
spiral for the 380-day mission. The M^/Mq 
can be seen to increase from 0. 467 to 0. 534 at 
a P/Mq of 150 watts per kilogram. The ac- 
tual change in mass ratio must take into ac- 
count the boost vehicle requirements of the no- 
initial- spiral trip, since for this case Mq is 
no longer the initial mass leaving the Earth orbit. 

Omission of the terminal spiral is accomplished by assuming that the returning vehi- 
cle can either be captured in a high energy ellipse or use atmospheric braking to the sur- 
face from an entry velocity equal to escape velocity (vk Circular “ 11 200 m/sec). 

For a mission analysis using this profile, provisions must be made for landing the men, 
such as including a reentry capsule in the payload or picking up the crew in the high en- 
ergy ellipse with a separate Earth- launched craft. Figure 5 shows the effects of omitting 
the final spiral for the 380-day mission. The M F /M Q increases from 0. 467 to 0 . 523 at 
150 watts per kilogram. Omission of the last spiral allows more time (of the order of 
20 days) for the heliocentric transfers, hence higher M^/Mq. 

Eliminations of both the first and last spirals could be brought about by a boost to 
escape energy at the start of the mission and atmospheric reentry from escape energy at 
Earth return. In figure 5, Mp/Mg is seen to increase from 0. 467 to 0. 600 at a P/Mq 
of 150 watts per kilogram, when both Earth spirals are omitted. Again, a boost vehicle 
must be taken into account and a reentry vehicle included in the payload. 


One-Thousand-Day Trip Time 

Figure C6 presents maps for 1000-day trips with wait times of 300, 400, and 500 
days. Figures 6(a) and (b) show how the final mass ratio increases with decreasing wait 
time. In figure 6(a) the maximum final mass had been recorded at each P/Mq by finding 
the best F/Mq. Along each dashed curve in figure 6(b), however, P/Mq and F/Mq are 
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Maximum final mass ratio, (Mp/M 0 ) max Final mass ratio, Mp/M 0 Maximum final mass ratio, (Mp/M 0 ) 


t 



Power to initial mass ratio, P/Mq, W/kg 

(a) Optimum thrust to initial mass ratio. 



Figure 6. - Effect of wait time for 1000-day 4-spiral trip. 



Figure 7. - Long against short mission times. Maximum 
final mass ratio; four-spiral missions. 


held constant at different arbitrary values. 

The solid curve, for contrast, has P/Mq 
fixed at 90 watts per kilogram, but the opti- 
mum (highest M f /M q ) value of F/Mq has 
been found at each wait time. Calculation dif- 
ficulties were too great to allow extension of 
the data below 200 days of wait for each P/Mq 
and F/Mq pair. 

Decreasing wait time for a fixed mission 
time forces the heliocentric portions of the 
trip to use higher elapsed times and travel an- 
gles. The need for return leg travel angles 
greater than 37r/2 radians contributes to the 
difficulties mentioned previously. The propel- 
lant consumption of low-thrust vehicles usually 
decreases with increasing heliocentric time. 

It is therefore probable for any fixed pair 
of P/Mq and F/Mq that M^/Mq will con- 
tinue to increase until wait time equals zero. 

The data presented in figures 6(a) and (b) 
are in direct contrast with the so-called im- 
pulsive Hohmann round trip. When impulsive 
trajectories are used with the stipulation that 
the velocity pulses must occur at the beginning 
and end of each trajectory, AV must increase 
if the travel time on a leg is found to be greater 
than the Hohmann time (260 days). Therefore, 
for a 1000-day mission, impulsive trips have 
an optimum wait time of about 450 days while 
it appears that low-thrust trips have an opti- 
mum of zero days. It might be advisable, how- 
ever, to stay at Mars for the longer wait times 
presented in the maps in order to do extensive 
investigation of the planet and avoid exposing 
the astronauts to the hazards of heliocentric 
transfer for an unnecessarily long time. 

Figure 7 compares (M f /Mq) for a 

1000-day trip with a 400-day wait and a 420- 
day trip with a 10-day wait. In general, longer 
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trips will allow the use of lower P/Mq, to result in the same final mass fraction, or 
will yield higher Mp/Mg values at the same P/Mg. Since powerplant fraction is propor- 
tional to P/Mg, increasing trip time appears, as expected, to benefit payload capability. 
For P/Mq = 90 watts per kilogram, the long trip time yields = 0. 735 as op- 
posed to 0. 441 for the short trip time. max 


CONCLUDING REMARKS 

Performance maps have been presented for constant low-thrust round-trip missions 
of short and long duration. For simplicity, it is assumed that no mass change occurred 
at Mars. It has been found that generation of these maps by the method used herein is 
very time consuming and would be far too costly to provide the large amount of trajectory 
data needed in extended mission analysis. The problem cases were those that (1) were 
near the all- propulsion boundary or (2) had travel angles on the return leg greater than 
37 t/ 2 radians, which is typical for trip times above 500 days with a waiting time of the 
order of 10 days and for 1000-day missions where the wait time is below 200 days. If, 
in the future, further improvements can be made in the methods of solving the variational 
trajectory problem for constant thrust, the computing time element may be of less con- 
cern. The maps are extremely useful as a basis for checking more approximate solutions 
to the trajectory problem. 

Since each point on a map represents an optimum round trip in terms of maximizing 
final mass ratio for a given P/Mq and F/Mg, there is no need to piece together one- 
way trips through a trial and error search. 

For constant- thrust trajectories, the all-propulsion boundary represents an impor- 
tant limitation on the area of valid solutions, that is, continuous thrusting on both legs of 
the heliocentric transfers. 

As opposed to high-thrust solutions, final mass ratio of low-thrust solutions in- 
creases monotonically with increasing mission time. 

A high-thrust, Hohmann-type round trip to Mars possesses an optimum waiting time 
(about 450 days) at a mission time of approximately 1000 days. It is interesting to note 
that this characteristic of high-thrust trajectories does not appear in the case of low 
thrust. For both the 1000-day missions and the short missions presented herein, final 
mass is seen to improve with decreasing wait time. 

Also shown was that a change in mission profile from the four- spiral trip may re- 
sult in improved vehicle performance. The four- spiral mission (as shown in appendix B) 
requires a low value of specific powerplant mass a to be competitive with nuclear or 
chemical rockets at short mission times. For high a;, other mission profiles such as 
boosting beyond escape energy and reentry from hyperbolic velocity were studied with 
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approximate techniques in reference 7 and were shown to improve the competitive 
position of electric propulsion. Reference 7, whose trajectory solutions were checked 
against accurate performance maps, contains a more complete study of the effect of mis- 
sion profile. This report does not contain accurate solutions to these other profiles be- 
cause of the prohibitive cost of obtaining them. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 24, 1965. 
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APPENDIX A 
SYMBOLS 


e 

eccentricity of vehicle tra- 
jectory 

R 

distance of vehicle from Sun, 
m 

F 

thrust, N 

T 

elapsed time, sec 

g 

gravitational constant of 

V 

velocity 


Earth, 9.80665, m/sec^ 

a 

specific powerplant mass, 

T SP 

specific impulse of thrustor, 


kg/W 


sec 

V 

overall engine efficiency 

K 

switching function (see ref. 3) 

V-p 

power efficiency of engine 

It 

semilatus rectum of vehicle 
trajectory, m 

^u 

fraction of mass leaving vehi- 
cle providing thrust 

M 

mass of vehicle, kg 

e 

true anomaly of vehicle trajec- 

m f 

final mass of electric vehicle, 


tory, radians 


kg 

X V X 2’ X 3’ 

Lagrangian multipliers in cal- 

m l 

useful payload mass of elec- 
tric vehicle, kg 

X 4’ X 5 

culus of variations solution 
(see ref. 3) 

M p 

mass of electric vehicle pro- 

<P 

elapsed travel angle, radians 


pellant, kg 

* 

thrust angle relative to local 

M PP 

mass of electric vehicle 


horizontal, radians 

powerplant, kg 


argument of pericenter, angle 

m tot 

total mass flow rate 


measured from starting 

M 0 

P 

P T 

initial mass of electric 
vehicle, kg 

effective jet power, W 

total power output of power- 
plant, W 


point of each trajectory to 
pericenter radius in the di- 
rection of motion, radians 
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APPENDIX B 


USE OF PERFORMANCE MAP FOR MISSION ANALYSIS 

The simplest concept of an electric vehicle is that it consists of three integral parts, 
propellant, powerplant, and payload. If this is the case and if a fixed payload and struc- 
ture are assumed, the difference between the initial mass Mq and the final mass Mp 
can be regarded as the propellant Mp necessary to fly the mission: 

M p = M q - M F (Bl) 

The powerplant mass Mpp is related to the total power output Pp by the specific 
powerplant mass a: 


Mpp — aPp 


(B2) 


The payload M^ is then the initial mass minus the propellant and powerplant 
masses: 


Mjj — Mq - Mp - Mpp 


or in terms of Mp and aPp, 


Mp Mp — QfPp 

Normalizing equation (B4) with respect to Mq yields 



(B3) 


(B4) 


(B5) 


In this simple analysis, when inefficiencies in the thrustor are ignored, Pp/M Q is 
equivalent to P/Mq appearing as a parameter on the performance maps. Figure 8 is a 
typical performance map showing one line of constant P/Mq. Since the performance 
map presents M^/Mq along lines of constant P/Mq, M l /M q can readily be obtained. 
This process is illustrated in figure 8 where a has been taken to be 3x10” ^ kilogram 
per watt. When thrustor inefficiencies are ignored, o' P/Mq is constant along a P/Mq 
line and the M l /Mq line is vertically displaced by a uniform amount below the M f /Mq 
line; hence, maximum M f /Mq gives maximum M l /Mq. 
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Maximum payload ratio. Mg/Mg Mass ratio, M/Mq 



Figure 8. - Final and payload mass ratios against thrust to initial 
mass ratio. Mission time, 380 days; wait time, 10 days; specific 
powerplant mass, 3. 0xl0~ 3 kilogram per watt; effective power to 
mass ratio, 90 watts per kilogram. 


Repeating this process for all P/Mq 
on the full map (fig. C2(b)) and then plot- 
ting maximum M^/Mq against P/Mq 
yield the overall maximum M^/Mq and 
optimum P/Mq. This is illustrated in 
figure 9. 

To determine the optimum thrust to 
mass ratio F/Mq, one plots the F/Mq 
that gives maximum MjVMq at each 
P/Mq against P/Mq and from the opti- 
mum P/Mq determines optimum F/Mq. 

Optimum specific impulse Igp is 
then found from the optimum values of 
P/Mq and F/Mq from the relation 



_P 

M 


0 


l SP 


JF 

M 


M, 


(B6) 


TOT £ 


0 


It should be noted that the optimized 
parameters are for a fixed value of a, 
and an entirely new set would be found 
by repeating the procedure just followed 
for a new a. 

Figure 10 shows optimum specific 
impulse, thrust to mass ratio, power to 
mass ratio, and maximum payload frac- 
tion as a function of a for a 380-day, 
10-day wait, four- spiral mission. The 

F/Mq, and P/Mq are all seen to decrease with increasing a, and M^/Mq rapidly 
approaches zero. 

Having a set of performance maps for different trip times but constant wait time is 
thus sufficient to determine payload ratio against trip time over a range of a. 


figure 9. - Maximum payload ratio at each power to mass ratio. 
Mission time, 380 days,- wait time, 10 days; specific power plant 
mass, 3. 0xl0"3 kilogram per watt. 
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The specific impulse of the thrustor Igp 
weight flow rate of ions contributing to thrust. If 7] 


is defined as the ratio of thrust to 
represents that fraction of mass 


u 


leaving the vehicle which is providing thrust (the propellant utilization efficiency), 
Igp ENG may be expressed as 
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Specific powerplant mass, a, kg/W 


Figure 10. - Optimum specific impulse, thrust to initial mass ratio, 
power to initial mass ratio, and maximum payload ratio against specific 
powerplant mass. Mission time, 380 days; wait time, 10 days. 
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For electric thrustors, the power efficiency of the engine rjp can be presented as 
in the sketch as a function of the specific impulse of the thrustor. 

If rj u is known, it is possible at each point on the 
map to determine the overall engine efficiency tj 
(n = 7 j u ?]p) and the total power requirements when effi- 
ciency is taken into account: 
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where rjp = ^pOsP, ENG>’ 


i n ii ii ia i 





Equation (B8) may not be of interest if a detailed study of the effect of the individual 
efficiencies r) u and 77 p is not required. A more simple analysis, showing the effect of 
overall efficiency, can be made when the variation of 77 is known or can be assumed for 
an electric thrustor system directly as a function of F/Mp^p. In this case, 


P T _ 1 / p\ 

M 0 * \ M 0/ 


(B 9 ) 


where 17 = 77(F/M TOT ) = i?(Ig p ), and Igp is given by equation (B6). 

Thus, at each point on the map the total power requirements are given, and the ef- 
fects of efficiency can be determined without the need for recomputing the trajectories 
since F was interpreted as total thrust and as all the mass leaving the vehicle. 

A further discussion of the effects of efficiency, structure, boosting, atmospheric 
braking, etc. , on payload fraction is beyond the scope of this report. References 7 and 9 , 
which are concerned with mission analysis, more fully treat this subject. The data of 
reference 9 were generated with a one-way calculus of variations program similar to the 
round-trip program employed here. 
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APPENDIX C 

PERFORMANCE MAPS 

This appendix presents more detailed data on the performance of the constant-thrust, 
constant-specific-impulse vehicle for the Mars round-trip mission. 

Table Cl presents a time history of the vehicle trajectory for a 380-day, 30-day-wait 
mission with no initial or final Earth spirals. 

Figure Cl(a) is a polar coordinate plot of this trajectory. Figure Cl(b) plots the tra- 
jectory for a 1000-day four- spiral trajectory with 400 days of wait. 

Figures C2 and C3 are performance maps for short mission time, four- spiral mis- 
sions, with various wait times at Mars. 

Figures C4 and C5 are maps for short mission times where the initial and/or ter- 
minal Earth spirals have been deleted. 

Figure C6 presents data for a 1000-day, four- spiral mission for three different wait- 
ing times at Mars. 
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TABLE Cl. - TRAJECTORY VARIABLES AS FUNCTION OF TIME FOR 380-DAY TRIP 


[Thrust to initial mass ratio, F/Mq, 2. 2555*10" J NAg; effective power to mass ratio, P/Mq, 100 W/kg; 
Lagrangian multiplier, X^, 267.93; constant of integration, 1.9229x10"^.] 


(a) Outbound leg 


Elapsed 

time, 

T, 

Elapsed 

travel 

angle, 

Distance of 
vehicle 
from Sun, 

Mass of 
vehicle, 
M, 

Thrust angle 
relative to 
local 

Switching 

function, 

K 

Lagrj 

ngian nniltiplic 
of variations 

rs in calculus 
solution 


Semilatus 
rectum of 
vehicle 

Eccen- 
tricity of 
vehicle 

True 
anomaly 
of vehicle 

Argument 
of peri- 
center, 

sec 

<P, 

radians 

R, 

m 

kg 

horizontal, 

radians 


X 1 

X 2 

X 3 

*5 

trajectory, 

f, 

m 

trajectory, 

e 

trajectory, 

Q, 

radians 

tu, 

radians 

0 

0 

1. 4950CX10 11 

1000. 0 

0.6386 

5.4623 

4.3441X10" 2 

5.8518X10" 3 

2. 5836X10" 8 

1.0000 

1.4950X10 11 

0 

0 

0 

1.7280X10 5 

. 03458 

1.4952 

995.60 

.6613 

5.2045 

4.2976 

5. 5227 

2. 5060 

1.0280 

1. 5264 

.02249 

. 3822 

-.3476 

6.9120 

. 14023 

1.4988 

982.42 

.7361 

4.4364 

4. 1232 

4. 5510 

2.2683 

1. 1065 

1.6198 

. 09044 

.4667 

-.3264 

1.3824X10 6 

. 28429 

1. 5121 

964.84 

.8525 

3. 4358 

3.7965 

3.3186 

1. 9509 

1. 1983 

1.7401 

. 18129 

. 5884 

-.3042 

2. 0736 

. 42947 

1. 5380 

947.25 

.9893 

2.4779 

3.3506 

2.2023 

1.6547 

1.2753 

1.8510 

. 27072 

.7202 

-. 2908 

2. 7648 

. 57248 

1. 5788 

929.67 

1. 1500 

1. 5736 

2.7861 

1.2467 

1.4048 

1.3377 

1.9453 

.35651 

.8618 

-.2893 

3.4560 

.70981 

1.6360 

912.09 

1.3458 

.7259 

2. 1172 

.4846 

1.2196 

1. 3855 

2.0124 

.43574 

1.0144 

-.3046 

4. 0903 
Coast 

.82795 

1.7033 

895. 95 

1. 5920 

0 

1.4314 

-.0304 

1. 1128 

1. 4169 

2.0349 

. 49868 

1. 1698 

3418 

6. 3833 

1. 1804 

1. 9868 

895. 95 

1. 5920 

1 o 

-1.3288 

-.5343 

1. 1038 

1. 4169 

2.0349 

. 49868 

1. 5222 

3418 

6.7789 

1. 2247 

2.0293 

887. 16 

-1.8236 

. 3828 

-1.7 593 

-.4545 

1. 1333 

1.4334 

1. 9926 

.46112 

1.6100 

-. 3853 

7. 4201 

1. 3072 

2. 1051 

869. 58 

-1.6455 

! 1.2027 

-2.6218 

-.1961 

1. 2035 

1.4783 

1. 9485 

. 39162 

1.7619 

-.4548 

8. 1113 

1.3840 

2. 1676 

852.00 

-1. 5178 

2.0889 

-3. 4844 

. 1848 

1.2868 

1 5426 

1.9462 

.32444 

1.8909 

-. 5070 

8.8025 

1.4572 

2. 2161 

834. 42 

-1.4150 

j 3.0444 

-4.3446 

.6825 

1.3836 

1.6291 

1. 9793 

. 25514 

2.0028 

-. 5456 

9.4937 

1. 5286 

2. 2504 

816. 83 

-1. 3265 

4.0757 

-5 1989 

1. 2959 

1.4957 

1.7406 

2.0451 

. 18037 

2. 1013 

57 26 

1.0185X10 7 

1. 5997 

2.2709 

799.25 

-1. 2469 

5. 1898 

-6.0418 

2.0281 

1.6268 

1.8807 

2. 1426 

.09742 

2. 1896 

-. 5899 

1.0717 

1.6551 

2.2775 

785.72 

-1. 1896 

6. 1089 

-6.6780 

2.6766 

1.7436 

2.0103 

2. 2394 

.02651 

2.2525 

- . 5974 

1.0904 
1. 1906 
Wait at Mars 

1.6748 

1.7809 

2.2779 

2.2779 

780. 97 
755.47 

-1.1700 

Spiral 

6.4444 

-6.8973 

2. 9223 

1.7885 

2.0607 

2. 2779 

0 

2. 3846 

-.7098 


(b) Return leg 


1. 4501X10 7 
1. 5469 

2.0556 
2. 1581 

2.2779X10 11 
2. 2779 

755.47 

730.84 

Spiral 
-2. 1244 

6.4444 

-6.0663X10" 2 

-3.7493X10 -3 

-1. 5067xl0" 8 

2. 2079 

2.2779X10 11 

0 

-3. 1416 

3. 1416 

1. 5642 

2. 1763 

2. 2775 

726. 45 

-2. 1094 

6. 1897 

-5.9421 

-3. 5506 

-1.4676 

2. 2595 

2. 2257 

.02960 

-2. 4469 

2.4651 

1.6160 

2.2298 

2.2713 

713.26 

-2.0640 

5.4373 

-5. 5591 

-2. 9879 

-1.3591 

2.4087 

2.0788 

. 11462 

-2.4033 

2. 4749 

1.6851 

2.2991 

2.2504 

695.68 

-2.0014 

4.4626 

-5.0314 

-2.3111 

-1. 2318 

2. 5946 

1.9057 

. 21946 

-2.3433 

2.4842 

1.7543 

2.3673 

2.2135 

678. 10 

-1. 9347 

3. 5217 

-4.4931 

-1.7113 

-1. 1200 

2.7653 

1.7598 

.31517 

-2.2790 

2.4882 

1.8234 

2.4358 

2. 1592 

660. 51 

-1.8615 

2.6164 

-3. 9514 

-1. 1823 

-1.0196 

2. 9206 

1.6422 

.40285 

-2. 2073 

2. 4850 

1. 96 16 

2. 5808 

1.9926 

625.35 

-1.6808 

. 9227 

-2.8784 

-.3179 

-.8379 

3. 1835 

1. 5004 

. 56221 

-2. 0257 

2. 4484 

2. 0307 

2.6627 

1.8782 

607.77 

-1. 5615 

. 1420 

-2.3529 

.0218 

-.7485 

3.2909 

1.4830 

.64096 

-1.9053 

2.4099 

2.0438 

2.6793 

1.8542 

604.44 

-1. 5358 

0 

-2.2545 

.0790 

-.7312 

3.3094 

1.4843 

.65643 

-1.8796 

2.4008 

Coast 

2.8482 

5. 5129 

1.0759 

604.44 

1. 5358 

0 

-1.3353 

-1.8182 

-.4699 

3. 3094 

1.4843 

.65643 

.9541 

2.4008 

2. 8654 

5.5771 

1. 1038 

600.05 

-2.4283 

.0368 

-1. 4924 

-1.7248 

-.3443 

3. 3338 

1.4477 

.62129 

1.0455 

2.3135 

2. 9000 

5.6938 

1. 1601 

591.26 

-2.2687 

. 1254 

-1.7933 

-1. 5039 

1082 

3. 3845 

1.3799 

. 55678 

1. 2236 

2.3122 

2. 9346 

5.7975 

1.2152 

582.47 

-2. 1102 

.2511 

-2.0794 

-1.2449 

. 1051 

3.4384 

1. 3216 

. 50005 

1.3948 

2.2446 

2.9691 

5.8905 

1.2673 

573.68 

-1. 9557 

.4308 

-2.3548 

-. 9538 

.2970 

3.4962 

1. 2750 

. 45040 

1. 5573 

2. 1750 

3.0037 

5.9750 

1.3152 
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-1.8081 

.6771 

-2.6231 

-.6345 

.4713 

3. 5591 

1.2418 

. 40622 

1.7086 

2. 1083 

3.0382 

6.0531 

1.3579 

556.09 

-1.6704 

. 9978 

-2.8870 

-.2886 

.6331 

3.6285 

1. 2230 
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1.8464 

2. 0487 

3.0901 

6. 1617 

1.4114 

542. 91 

-1.4859 

1.6250 

-3.7780 

.2790 

.8627 

3.7482 

1.2222 

.30493 

2.0758 

1. 9778 

3. 1592 

6.2972 

1.4617 

525.32 

-1.2809 

2.7326 

-3.7914 

1. 1310 

1. 1678 

3.9457 

1.2714 

.21572 

2.2187 

1. 9204 

3.2283 

6.4296 

1.4886 

507.74 

-1. 1145 

4. 1374 

-4.2859 

2. 1038 

1. 5014 

4. 2004 

1.3749 

. 10637 

2.3721 

1.8994 

3.2733 

6. 5172 

1.4946 

496.31 

-1.0211 

5. 2038 

-4. 5856 

2. 8098 

1.7485 

4.4046 

1. 4702 

.02107 

2.4553 

1. 9037 

3.2837 

6. 5318 

1.4948 

493.66 

-1.0007 

5.4671 , 

-4.6511 

2.9820 

1.8102 

4.4568 

1.4954 

0 

0 
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(a) Mission time, 380 days; wait time at Mars, 30 days; no initial or terminal Earth spirals; power to initial mass ratio, 
100 watts per kilogram; thrust to initial mass ratio, 2. 2555x10“-* newton per kilogram. 

Figure Cl. - Distance of spacecraft from Sun against angle traveled. 
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270 


(b) Mission time, 1000 days; wait time at Mars, 400 days; 4-spiral mission; power to initial mass ratio, 55. 78 watts per kilo- 
gram; thrust to initial mass ratio, 1. 185x10"^ newton per kilogram. 

Figure Cl. - Concluded. 
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(b) No initial spiral. 

Final mass ratio against thrust to initial mass ratio. Mission time. 340 days; wait time, 30 days 





Thrust to initial mass ratio, F/Mq, N/kg 
(d) No initial or final spirals. 
Figure C4. - Concluded. 





Final mass ratio, Mp/M 0 



(a) Four spirals. 



Thrust to initial mass ratio, F/Mq, N/kg 
(b) No initial spiral. 

Figure C5. - Final mission ratio against thrust to initial mass ratio. Mission time, 380 days; wait time, 30 days. 
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